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I. INTRODUCTION 

The concept of the brane and brane world exhibits broad possibihties in the particle 
theory beyond the standard model. In the early stage of study on the higher dimensional 
theories, it has been considered that we may live in a domain wall perpendicular to the 
fifth dimension [1] . The domain wall may be described by a kink configuration of an self- 
interacting scalar field. The matter field which couples to the kink is localized in the extra 
dimension and acquires a mass spectrum consisting both discrete and continuum modes in 
general. 

Even if we consider the model with a compact dimension, the kink or solitonic configura- 
tion is expected to play important roles. The localization of fermions due to solitonic domain 
walls can explain the hierarchy in masses and the Yukawa couplings [2]. The compactified 
dimension is taken as S^ or S^ /Z2 in such a case. The kink solutions in the periodic compact 
space are known to be topologically unstable like sphalerons. They have been studied by 
some authors [3]. 

In the recent decade, another approach to the particle physics models has been explored. 
The idea of dimensional deconstruction [4] is to construct the models mimicking the higher 
dimensional gauge theory from the four dimensional gauge theory. In some sense, the de- 
construction means introducing the discrete extra space into the theory. Some merits of 
higher dimensional theory can be inherited by four dimensional theory, whose high-energy 
behavior can be well controlled. 

Therefore it is worth examining the 'deconstructed kink model' in the framework of 
four dimensional theory. Fortunately, discrete models for kink-like solutions have been 
considered apart from particle physics [5]. Thus we must study the solutions with specific 
boundary conditions and especially coupling to fermions to apply the discrete models to 
particle physics. 

Most general starting point is to consider the multi-scalar model in four dimensions, 
whose Lagrangian density reads 

1 ^ 

^=^Y. ^'<t>^^.ct>^ - ^(01, 02, . . • , M , (1.1) 

^ 1=1 

where (j)i are real scalar fields and V is the potential. Unlike the original dimensional 
deconstruction scheme, the gauge invariance is not assumed. Moreover even if the continuum 



limit is uniquely defined, it is known that the discretized action can be exhibited by various 
different forms [5]. Since the other guiding principles are needed for detailed study, we will 
refer the supersymmetric model as shown later. 

The original dimensional deconstruction is found to be associated with the graph struc- 
ture [6]. Thus we begin with considering the bilinear term in scalar fields and construct 
the interaction terms on vertices or edges of a graph taking the correspondence with the 
continuum theory into consideration. 

In the next section, we review the kink solution in (1 + l)-dimensional field theory. 
Since we must include fermions later, we consider a supersymmetric theory; the structure 
of this theory will be referred in the later sections. In Sec. Ill, we consider discretized or 
deconstructed free scalar models, whose continuum limit is the compactified theory with 
S^ or S^/Z2- In Sec. IV, we choose the form of the scalar interactions, which is expected 
to make the 'discrete kink' configurations possible. The result of numerical calculations on 
the solitonic solutions is shown in Sec. V. In Sec. VI, we consider the Yukawa coupling 
of fermions and scalars. Various 'discrete' versions of the Yukawa coupling are possible, 
though their 'continuum' limits are all the same one. In this paper, we consider the models 
whose couplings are similar to those of the supersymmetric model mentioned in Sec. II in 
the continuum limit. The numerical result for the fermionic spectra is exhibited in Sec. VII. 
The last section is devoted to summary and prospects. 

II. KINK SOLUTIONS IN THE CONTINUUM THEORY 

In this section, we review a continuum theory for kink solutions [7]. The form of the 
action is used to determine the way for discretization of bosonic as well as fermionic sector 
later. 

We consider the (1 + l)-dimensional supersymmetric field theory, whose Lagrangian den- 
sity is 

^ = 2<^^' - 2 [-1^) + 2'^'W - 5^V* ■ (2-1) 

where we choose the superpotential W{(j)) as 

Wi<i>) = 9{j- «'</') , (2.2) 

with g and a are constants. 



Then the bosonic part of the Lagrangian density reads 

CB = \m'-\9'{<t>'-a')\ (2.3) 

and thus we find the equation of motion as 

-dj4> + dl4> - 2/(02 - a2)0 = . (2.4) 

A non trivial static solution, a kink solution, satisfies the equation of motion, and ex- 
pressed as 

0o(^) = atanh^'ax , (2.5) 

where we set an integration constant for (/>o(0) = 0. One can find that this solution also 
satisfies 

dA + 9{<f'-a^)=dA + ^ = ^. (2.6) 

The coupling of the scalar and the (Majorana) fermion is proportional to —g(j)il)ip. In 
the present case, it is not surprising and is of the usual Yukawa form. This coupling with 
general coupling constant (say, —Gcpipilj) or similar type of models have been investigated 
in the soliton physics [8] . 

The solitonic configurations in the periodic compact space have also been considered in 
several authors [3] . It has been shown that he nature of quasi-stable solutions depends on 
the size of the compact space. A solution of (2.6) is not necessarily a solution of the equation 
of motion in the compact space. Interestingly, these features are also true in some sense for 
the case of discrete solitons studied in Sec. IV. 

III. DECONSTRUCTING A FREE SCALAR FIELD THEORY AND A GRAPH 

Suppose A^ real scalar fields 0^ (i = 1, . . . , A^). The mass spectrum of the scalars is the 
same as that of the vector bosons in the dimensional deconstruction model [4] , if the action 
is given by 

i=l i 

N -, -, 

= Y.-oi^<t>^f - -J'Y.(<t>^+^- <l>^f ^ (3-1) 

1=1 i 



where A(j)i = 0^+1 — (j)i and / is a mass scale. In the four-dimensional point of view, we can 
say that the potential is of a 'difference-squared' type. 

The sum in the potential is interpreted as (J)n^i = 0i and from i = 1 to N . Then the 
large A^ limit leaving L = N / f constant leads to the five dimensional scalar theory with S^ 
(Kaluza-Klein) compactification. Namely, the N x N mass-square matrix in this model is 
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= fA{C) , 



(3.2) 



and the eigenvalues are 



4/^sin2^ (fc = 0,l,...,iV-l) 



(3.3) 



which becomes (27r/c/L)^ in the large A^ limit mentioned above. Here L corresponds to the 
circumference of S^. 

The other summation rule is possible. The sum in the expression (3.1) is taken from 
z = 1 to A^ — 1. In this case the N x N mass-square matrix becomes 
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rA(p) 



(3.4) 



and its eigenvalues are 



4/2 sin^ 



2N 



(fc = 0,l,...,iV-l) 



(3.5) 



and the large A^ limit yields the theory with the compactification on S^ /Z2. 

In conclusion so far, the construction of the potential is associated with the structure of 
a graph, in the sense of graph theory [9]. A graph consists of vertices and edges each of 
which connects two vertices. The edge e connects the vertices o(e) and t(e). The vertex o(e) 
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means the origin of the edge e while the vertex t(e) means the terminus of the edge e. We 
assign A^ scalar fields onto A^ vertices of the graph and their action is written by 



where V and £ denote the set of vertices and edges, respectively. 

Now we can express the mass term by using the incidence matrix E defined as 



(3.6) 



iE\ 



1 if 'U = o(e) 
-1 if^ = t(e) , 
otherwise 



(3.7) 



since in both cases, the difference Acj) is defined as associated with each edge such that 



and then 
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(3.8) 



(3.9) 



where the matrix E^ is the transposed matrix of E. 

For the first case introduced in this section, the incidence matrix is 
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(3.10) 



and its transposed matrix is 



'^ 1 -1 
1 -1 



E^iO 



•■ \ 
•■0 



(3.11) 



■•■ 1 -1 

\^-i ■•■ 1 y 

The corresponding graph is the cycle graph in the term of graph theory (Fig. 1). Note that 




FIG. 1. The cycle graph C4 and the path graph P4. Each suffix indicates the number of vertices 
in the graph. 

A((7) = E{C)E^{C)^ where A((7) is known as a graph Laplacian for a cycle graph. 
For the second case introduced in this section, the incidence matrix is 
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(3.12) 



and its transposed matrix is 



E^{P) = 
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(3.13) 



\ w w w • • • 1 —1 / 

The corresponding graph is the path graph in the term of graph theory (Fig. 1). Note that 
A(P) = E{P)E^{P)^ where A(P) is known as a graph Laplacian for a path graph. 

In the next section, we will examine the nonlinear interaction of scalars, which is essential 
for the kink-like configurations, by referring the graph structure. 



IV. DECONSTRUCTING THE SCALAR POTENTIAL 



In the previous section, the term corresponding to the kinetic term of the scalar theory 
has been simply discretized. However, we should consider some constraints to the interaction 



term in the discretized action, because there are many possible descretizations which have 
the same continuum hmit. 

First of aU, we want to avoid apparently complicated interactions. We postulate renor- 
malizability, at least superficially. At most four scalar fields coupled mutually at a point in 
four dimensions. Secondly, we wish to consider models possessing a continuum limit for a 
large number of four-dimensional scalar field. In this limit, the static configuration is much 
alike kinks or sphalerons on S^ /Z2 or S^. Thirdly, the structure of the interaction should 
reflect a graph structure. The deconstructed 'kinetic term' is expressed as a sum over edges 
of a graph as shown in the previous section, so the deconstructed interaction term should 
also be written as a sum of the contribution assigned at edges. 

Now the possible Lagrangian density is the following: 






^t{e) + Q^t{e)(t^o{e) + P^t^lje) 2^ 



2p + q J 

(4.1) 

where p and q are dimensionless constants and we consider only a path graph or a cycle 

graph, corresponding to the continuum theory on S^ /Z2 or S'^, respectively. 

For a cycle or path graph, the potential term in the sense of the four-dimensional theory 

is 

Setting the derivative with respect to ^^^^^ {k is not 1 or A^ for the case with a path graph) 
to zero yields the recurrence relation among three terms. Unfortunately, it seems that this 
recurrence relation cannot be satisfied by multiple use of some recurrence relations in two 
terms. Thus for a finite number of fields, similar relations to the continuum theory like (2.6) 
cannot hold. 

Because this is still far from a specific model due to arbitrary numbers p and g, we select 
the case by special conditions p = q = 1. This selection was adopted first by Speight and 
Ward (in [5]). In this case, we find 

= 9 (^^^ - "H'l'.+i - *)) = W''(*+i) - »'(*) ■ (4.3) 

This relation shows that the derivative of the superpotential ^ is replaced by the difference 



AW 

A<j> 



, or 



AW = -E^W , where W 



' H/(0i) ^ 



Thus in this case, we can say that the interaction terms is related to the graph structure. 

For p = q = 1, which we concentrate on this case hereafter, the potential minimum is 
given by the simultaneous equations 



-K— = -J {(Pk+1 - 20fc + (pk-ij + g I a I 



(4.4) 



3 / 3 

where k is not 1 or A^ for the case with a path graph. To obtain static scalar configurations, 
we should solve the recursion relation (4.4) for a cycle graph, and for a path graph with 
additional 'boundary' equations. 



dV „o, , ,. . 2 ( 4>l + 4>2(l)i + 4>\ 2\4>2 + 2(i) 



-r{h-<t^i) + 9' 



i = o, 



dV 



d(t> 



N 



-/^(-07v + 07v-i) + g 



2 f^N + 4^N(pN-l + (p%-i 



— a 



^N + 07V- 1 



(4.5) 



0. (4.6) 



NUMERICAL SOLUTIONS 



We can look for symmetric configurations by numerical calculations. The recursion rela- 
tion (4.4) can be rewritten as 



Vk+i - "^Vk + Vk-i , (yl^i + Vk+iVk + yl , ^ Vk+i + 2?/^ 



/l2 

^yl + VkVk-i + yl-i 



1 



+ 



^ I 2j/fe + Vk-i ^ Q 



(5.1) 



3 7 3 

where yi = (j)i/a and h = ga/ f . 

Although the third order equations can be analytically solved, here we use Mathematica 
and the command FindMimimum to solve the solution for potential minima, without detailed 
analysis on numerical errors. This is partly because our model is still a toy model and 
because we consider only modest number of fields as a considerable field theory. 

We try to find non trivial solution for A^ = 4, 8, and 16 for both cases with path and cycle 
graphs. 
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A. a soliton on a path graph 



FIG. 2 shows the sohtonic profile for A^ = 4, 8, and 16. Even for the smaU number of A^, 
a non trivial configuration exists for a sufficiently large value for h. Larger values of h seem 
to induce the excess of the value of 0^ beyond a, near the center of the kink. 



FIG. 2. Numerical solutions for the scalar field configuration in the cases with P4 for the left 
column, Ps foi" the center column, and Piq for the right column. Each vertical axis indicates 
Hi — 4^i/ci- The figures in the five rows correspond to the parameter h = 2"-*^, 2~^-^, 1, 2^-^, 2-*^, from 
the top to the bottom. 



B. a soliton on a cycle graph 



FIG. 3 shows the solitonic profile for A^ = 4, 8, and 16. The feature is the same as the 
case with the path graph. 

In both cases, a large h means that the nonlinear part of the potential dominates the 
difference-square part. The tendency indicates that the simultaneous limit A^ ^ oc and 
h ^ leads to the continuous solitons. 
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FIG. 3. Numerical solutions for the scalar field configuration in the cases with C4 for the left 
column, Cs for the center column, and Ciq for the right column. Each vertical axis indicates 
Hi — 4^i/ci- The figures in the five rows correspond to the parameter h = 2~^, 2~^-^, 1, 2^-^, 2^, from 
the top to the bottom. 

VI. YUKAWA COUPLINGS 



Now we consider fermion fields coupled to the scalar fields. Our starting point is the 
discretization of the kinetic term in the continuum action in Sec. II. It is already known 
[6] that the fermion kinetic term can be expressed in the incidence matrix (E). One set of 
chiral fermion fields ip are assigned on vertices of a graph and the other set of chiral fermion 
fields X with the opposite chirality are assigned on edges of a graph. Then the discretized 
kinetic term has the form 

-f^Ex + h.c. , 

where indices are suppressed. 

Now we postulate three possibilities of Yukawa interaction. We call them Case A, Case 
B, and Case C, for later convenience. 
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A. a definition from difference operations 

In the continuum theory, we have already seen that the scalar (massy matrix can be 
written in the form 

(6.1) 



del) 
in the two-dimensional (j)"^ theory in Sec. II, while in our model, 

AW ^ 



A^ 



gives the potential, where 



W^(0i) = '^-a\ 



In the supersymmetric model, the fermion mass matrix is often given in the form 



(6.2) 



(6.3) 



(6.4) 



in the continuum theory. The naive discretization of this form will be formally expressed as 



A<f) 



2 ■ 



(6.5) 



In actual, since AW/A(f> is 



9 \ o ^ 



(6.6) 



the difference of two sequential terms becomes 



i+l "T <Vi+l<i 



'i-l 
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•H+1 -r ^>i 



9- 



3 ''3 

k+i -<Pi + <Pi- 4>i-i) = 9- 



i+i + 0i + (^i_i , , , , \ 4>i+i + 4>i + 4>i-i 



+1 - <^i_i) , 

(Zi0i + Z\(/)i_i) . (6.7) 



Graph theoretically speaking, the difference A(\) lives on each edge of a graph. Therefore we 
force chiral fermion fields on edges to couple the second difference in (6.7). Namely, for a 
cycle graph. 



y^c) = ^ 
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(6. 
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couples to X as Yd{C)x, where 



X 



X2 



(6.9) 



\Xn J 
For a path graph, the couphng at both ends should be fixed by hand. We take 



YdiP) = o 



/ 201+02 

01 + 02 + 03 01+02+03 

02+03+04 02 + 03 + 04 



V 



\ 



/'7V-2 + 07V-1 + 07V 07V-2 + 07V-1 + 07V . 
07V-1 + 207V / 



which couples to x ^s y^(P)x, where 



X 



X2 



V Xn-\ j 
We consider the four-dimensional fermion bilinear term 



(6.10) 



(6.11) 



f^{-E^h'a-^Ya)x + h.c., 



(6.12) 



where ip = ('0i, '02, • ' ' ? V^a^) is fermions on vertices of the opposite chirahty to x- The 
couphng h' is Ga/ f ^ where G is the general Yukawa coupling. 

B. a definition from SUSY inspired structures 

Suppose a cycle graph. The superpotential Wi can be defined on each edge as 



wm.m) = m^ 



/([(/)],+! -[0]i)+^ 



'M?+i + Mm[0]i + 



— a 



(6.13) 



where the superfields represent multiplets, each of which consists of a scalar, a chiral fermion, 
and an auxiliary field: 

[0]i9(^i,t/;,,F/'), (6.14) 



[^,3{^,,Xi,Ft). 



(6.15) 
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Then we get 



'dm 



f(cPi+,-<P,y + g 






^2/,f^m-^? .. 



V 3 
Assuming $i = 0, the potential for 0j is given by 



- a [(pi+i - (pi) 



(6.16) 






2 AT 

I =E 

i=l 
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'^i+1 - (pi 



r+ 



>f+l + 0i+10i + 02 



— a 



(6.17) 



since there is no terminal vertex in a cycle graph. This is just the scalar sector in our model 
considered in the present paper. 

Thus we can derive the coupling terms from (6.13). The fermion bilinear term is found 
to be 



i 

Now we define 



f{^i+i-tpi)+g- 



2<j)i + <j)2 (/'i + 2(/.2 

2(f)2 + (f>3 (f>2 + 2(/>3 

1 



+ h.c. . 



(6.18) 





ycpN + 201 

and rewrite the biUnear term as 
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• • 20iv. 


-1 +0iV 
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1 + 20iv 










20 


W +01 



(6.19) 



fxi-E' + ha-%' )^ + h.c. . 

The couphng h may be changed to h' in a general case. 
In the case with a path graph, 

/ 



(6.20) 



r^(P) - 1 



201+02 01+202 

202 + 03 02 + 203 








(6.21) 



\ • • • 20Ar_i + 0Ar 0Ar_i + 20Ar J 

can be selected. Note in this case, although the terminal vertices may break an exact 
supersymmetry. Of course, supersymmetry is merely a guiding principle here. 
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C. a definition from a simple democracy on edges and vertices 

Another simple choice respecting a graph structure for the Yukawa couphngs is possible. 
The definition is 
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• (pN-1 4>N 
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■ (j)N 
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for a cycle graph and also 



K^iP) 



01 02 ••• 

02 03 ••• 







5JV-1 (pN I 



for a path graph. 

The four dimensional fermion bihnear term is 



(6.22) 



(6.23) 



(6.24) 



We considered three cases for Yukawa couplings. In the next section, we will obtain the 
spectrum of the fermions for each case. 

VII. MASS SPECTRA FOR FERMIONS AND LOCALIZATION 



FIG. 4 shows the fermion mass spectrum for the model associated with a path graph 
and that with a cycle graph. The eigenvalues of the matrix {—E + ha~^Y){—E^ + ha~^Y^) 
are plotted, where the sufiices indicating definitions of the Yukawa coupling are suppressed. 
The Yukawa coupling constant is taken as h^ = /i, so that the model in the continuum limit 
is expected to correspond with the supersymmetric model in Sec. II. 

For the mode with a path graph Pie, we can see the degeneracy of eigenvalues at /i ^ 1. 
In other cases, the degeneracy may occur at slightly large value for h. 

It is known that the fermion spectrum has the finite discrete modes and continuous modes 
in the background of a kink in the continuum theory [8] . The degeneracy of our case refiects 
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FIG. 4. The spectra of the fermion mass-square in the case with a path graph for the left column and 
a cycle graph for the rught column. The figures in the three rows correspond to A^ = 4, 8, 16, from 
the top to the bottom. The horizontal ticks correspond to the parameter h = 2~^, 2~^-^, 1, 2^-^, 2^, 
from the left to the right. The red points indicate the eigenvalues in Case A, the green points 
indicate the eigenvalues in Case B, while the blue points indicate the eigenvalues in Case C. 



the continuum spectrum in the corresponding field theory. This phenomenon bring about 
for the appropriate coupling /i, because the scalar configuration leaves far from the kink 
shape for small and large h. For a small /i, the difference-squared part is dominant. For a 
large /i, the non-linear potential is dominant. In both cases the shape of a kink is deformed. 

Finally we show the localization of zero mode fields. In applications to the particle 
physics, the localization of matter fields is expected and is substantial for inclusion of other 
interactions among the matter fields [1, 2]. In FIG. 5, the eigenvectors associated with zero 
eigenvalues for Piq with h = 1 are exhibited. 

The difference of profiles of zero modes is minute for the three choices for the Yukawa cou- 
plings. As in the model with a continuous extra dimension, we can utilize the 'localization' 
for constructing elaborated gauge interacting models with several gauge groups. 
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FIG. 5. The profiles of the eigenvectors associated with zero modes for Piq with /i = 1. The red 
points indicate the eigenmodes in Case A, the green points indicate the eigenmodes in Case B, 
while the blue points indicate the eigenmodes in Case C. 

VIII. SUMMARY AND OUTLOOK 

We investigated the solitonic configuration of scalar fields in deconstructed theory, or 
in a discrete extra space. It is expected that a large number of fields realizes the thee 
configuration corresponding to the kink in the continuum space. 

We found discrete eigenvalues for fermions appropriately coupled to the kink-like config- 
uration of scalar fields and the appearance of degeneracy of them for a set of certain critical 
couplings. This degeneracy is not observed in the continuum theory, where the continuum 
spectrum can be found. We also emphasis on the simple features of the models: the lo- 
calization of matter fields occurs and the existence of zero modes, which correspond to the 
standard-model particles, is guaranteed. 

The application of the deconstructed scalar model to particle physics turns to be in- 
teresting. The degeneracy in the spectrum may lead to the significant quantum eflFects 
on the zero mode fields if the other interactions are turned on. The diflFerence from the 
higher-dimensional model with compactification is expected to be found, since the contin- 
uum theory includes infinite tower of fields. Moreover the study of the Casimir eflFect will 
be important especially if gravitation is incorporated into the models. We think that the 
dynamical creation and evolution of the kink or 'topological' (not in the strict meanings) 
configuration are worth studying in the cosmological context. 

The generalization of models to those with the structure of general graphs other than cycle 
and path graphs is mathematically interesting. We wish also to investigate deconstructing 
non-linear sigma models and 'topological' configurations in them in the future works. 
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